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Abstract 

We study the asymptotic expansion for the Landau constants G n , 

OO 

7 T G n ~ ln(16AT) + 7 + ^ -^r as n oo, 
k =1 1 


where N = n + 1, and 7 is Euler’s constant. We show that the signs of the coefficients 

, \ 1(1 +1) 

ctk demonstrate a periodic behavior such that (—1) 2 au +1 < 0 for all l. We further 
prove a conjecture of Granath which states that (—1) £i(N) < 0 for l = 0,1,2, • • ■ and 

n = 0,1,2, ■ ■ •, £i(N) being the error due to truncation at the Lth order term. Consequently, 
we also obtain the sharp bounds up to arbitrary orders of the form 


ln(16iV) +y + J2w <7rGn < ln ( 16iV ) + 7 + X Tfk 


Ofc 


fc=1 


N k 


fc=1 


N k 


for all n = 0,1,2 • ■ •, all p = 4s + 1, 4s + 2 and q = 4m, 4m + 3, with s = 0,1, 2, • • • and 
to = 0 , 1, 2 , ■ ■ ■ . 
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1 Introduction and Statement of Results 

In 1913, Landau m proved that if f(z) is analytic in the unit disc, and \f(z)\ < 1 for \z\ < 1, 
with the Maclaurin expansion 

f(z ) = a 0 + a\z + a 2 z 2 -|-b a„z" -|-, \z\ < 1, 

then there exist constants G n such that 

|a 0 + ai + --- + a n | < G n , n = 0,1, 2, ■ ■ ■ , 
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National Natural Science Foundation of China under grant numbers 10871212 and 11571375. 
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and the bound is optimal for each n, where Go = 1 , and 


G 


n 



fl -3- ••• (2n- 1)\ 
V 2 • 4 • • • • ( 2 n) ) 


for n = 1 , 2 , • • • . 


( 1 . 1 ) 


The constants G n are termed the Landau constants. The large-n behavior is known from 
the very beginning. Landau m derived that 

G n ~ — In n as n —> oo; 

7r 


see also Watson [18] . It is worth mentioning that there exist generating functions for these 
constants; cf. [7], possible q -versions of the constants; cf. m , and an observation made by 
Ramannujan (cf. [7]) that relates the Landau constants to the generalized hypergeometric func¬ 
tions. Useful integral representations for G n have been obtained from such relations; cf., e.g., 
Watson p~8] : see also Cvijovic and Srivastava [7]. 

The approximation of G n has gone in two related directions. One is to obtain large-n 
asymptotic approximations for the constants, in a time period spanning from the early twentieth 
century Effing to very recently E3H23- The other direction is to find sharper bounds of G n for 
all nonnegative integers n. Authors working on the sharper bounds includes Brutman [3] and 
Falaleev [§] (in terms of elementary functions), Alzer [2] and Cvijovic and Klinowski [ 6 ] (using 
the digamma function), Zhao [20] , Mortici [IT] and Granath [9] (involving higher order terms), 
and Chen and Choi [5] and Chen [4] (digamma function and higher order terms). The list is by 
no means complete. The reader is referred to [[71 1121 fT3lj for a historic account. 


1.1 Optimal bounds up to all orders 

Attempts have been made to seek bounds in a sense optimal, and up to arbitrary accuracy. 

In 2012, Nemes [15] derived full asymptotic expansions. For 0 < h < 3/2, he shows that the 
Landau constants G n have the asymptotic expansion 

G n ^ — ln(n + h) + —(7 + 4In 2 ) — 7^~rT\k as n +°°’ ( 1 - 2 ) 

7 r 7 r (n + hr 

k> 1 v ' 

where 7 = 0.577215 • • • is Euler’s constant. Earlier in 2011, the special cases h = ^ and h = 1 
were established by Nemes and Nemes m using a formula in [16]. They also conjecture in m a 
symmetry property of the computable constant coefficients such that gk(h) = (— l) fc <7fc (3/2 — h) 
for every k > 1. The conjecture has been proved by G. Nemes himself in m ■ A natural 
consequence is that for h = 3/4, all odd terms in the expansion vanish. In this important 
special case, Nemes [15j has further proved that 

Proposition 1. (Nemes) The following asymptotic approximation holds: 

7 tG 

where the coefficients (—l) s_l_1 /32s are positive rational numbers. 

The derivation of Nemes m is based on an integral representation of G n involving a Gauss 
hypergeometric function in the integrand. An entirely different difference equation approach is 
applied in Li et al. [12] to obtain full asymptotic expansions with coefficients iteratively given. 


ln(n + 3/4) + 7 + 4 In 2 + ^ 


/?2 s 


7 ^ (n + 3/4) 


2s : 


n 


00 , 


(1.3) 
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What is more, in a follow-up paper M, it is shown that the error due to truncation of m 
is bounded in absolutely value by, and of the same sign as, the first neglected term for all 
n = 0,1, 2, • • • . An immediate corollary is 

Proposition 2. (Li, Liu, Xu and Zhao) For N = n + 3/4, it holds 

2m q 2k—1 q 

l niV + 7 + 41n2 + ^^< 7 rG n <lnfV + 7 + 41n2 + ]T ^ (1.4) 

S =1 S=1 

for all n = 0,1,2,• • • , m = 0,1,2, • • • , and k = 1,2, • • • . 

In a sense, the formulas m and (11.41) in the above propositions seem to have ended a 
journey since one has thus obtained optimal bounds up to arbitrary orders. Yet there is an 
interesting observation worth mentioning, as presented in the 2012 paper [9] of Granath; see 
also [13]. 

Granath derives an asymptotic expansion 

OO 

7 rG n ~ ln(161V )+7 + ^-^, n-)■ oo, (1.5) 

k =l 1 

where a*, are effectively computable constants but not explicitly given, except for the first few. 
Here and hereafter we use the notation N = n + 1. 

Denoting the truncation of (11.51) 

i 

A l (N) = ln(16N)+ 1 + J2^, ( 1 - 6 ) 

k =1 1 

then one of the main results in Zhao [2D] reads A 2 (N) < ixG n < A 3 (IV) for n = 0 , 1 , 2 , ••• . 
Mortici m have actually proved that A$(N) < nG n < A^(N) for all non-negative n. 

In [9], Granath proves that A^(N) < n G n < Aj(N) and states that Ag(N) < irG n < An(N), 
for all non-negative n. Based on these formulas and numerical evidences, Granath proposes a 
conjecture. 

Conjecture 1. (Granath) It holds 

i(i+ 1 ) 

(-1)— (nG n -Ai(N)) <0 (1.7) 

for all n = 0 ,1, 2 , • • • and l = 0 ,1, 2 , • • • . 


1.2 Statement of results 

We will show that the conjecture is true. To do so, we will make use of the second order 
difference equation for G n employed in and some estimating techniques used in [Hi- 
First we denote the error term 

ei(N) = 7 rG n - A(N) = v vG n - |ln(16AT) + 7 + J2 w] ’’ (L8) 

cf. (11.61) . where N = n + 1. It is readily seen that £i(N ) ~ ai + \/N l+1 as N —>■ 00 . Hence we 
may start by showing that (11.71) holds for large n. To this aim, we have 
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Theorem 1 . The coefficients of the asymptotic expansion (USD satisfy 

i(i+i) 

(- 1 ) — a l+1 <0, 1 = 0 , 1 , 2 ,--- . 

Next, we will prove the conjecture for all non-negative n. 

Theorem 2. For N = n + 1, it holds 
yi+ 1 ) 

(-1) — £l (N)<0 


(1.9) 


( 1 . 10 ) 


for l = 0 , 1 , 2 , • • • and n = 0 , 1 , 2 , • • • . 

As a straightforward application of Theorem ( 2 j we obtain the following sharp bounds up to 
arbitrary orders. 

Corollary 1 . For N = n + 1, it holds A p (N) < irG n < A q (N), that is, 

ln(161V) + 7 + fffik < nGn < ln ( 16A 0 + ( L11 ) 


for all n = 0,1, 2, • • • and for all p = 4s + 1, 4s + 2 and q = 4m, 4m + 3, with s = 0,1, 2, • • • 
and m = 0 , 1 , 2 , • • • . 


In view of Theorem [TJ we see that the bounds in (11.111) are optimal as n —>• oo. 

Theorem [2] can actually be understood as an estimate of the error term, such that the error 
due to truncation is bounded in absolute value by, and of the same sign as, the first one or two 
neglected terms. Indeed, since 

e,(N) = %&+e l+1 (JV) and e,(JV) = $& + $g + e,«W, 
taking into account the signs in Theorems Q] and El we have 


0 < S4k+i(N) < 


Q4/C+2 «4fc+3 

j\r4fc+2 _/\r4fc+3 


and 


0 < £4fc+2(A0 < 


«4fc+3 

jyik+3 


for all non-negative integers n and k, and 


Q4fc+1 

_/y4fc+i 


< £4 k{N) < 0 


and 


Q4fc+4 «4fc+5 
_/y4fc+4 + jy4fc+5 


< £4fc+3(A0 < 0 


for all non-negative integers n and fc. 

As a by-product of the proof of Theorem El we have approximations of the asymptotic 
coefficients, follows respectively from (|4.10l) and (|4.15l) : 

Corollary 2. Assume that a & are the coefficients in the asymptotic expansion (11.51) . Then 
we have 


, ,, l+ i 2 ( 2 fc- 2 )! 

“* = (_1) “(Spw 


1+0| r 


and 


, i \k+i 8 ( 2 fc)! ln( 2 fc + 1 ) Jl)) 

“ 2fc+1 = ( ” 1} (27r) 2fc + 2 V 1 + ° UM J 


( 1 . 12 ) 


(1.13) 


as k —>• oo. 


□ 
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2 The asymptotic coefficients and the proof of Theorem [T| 

From the representation (11.111 one obtains the recurrence relation 


Gn +1 G n — 


2 n+ 1 


(Gn - Gn- 1) . 


2 n + 2 / 

Set N = n + 1, we may rewrite it as a standard second-order difference equation 

1 


w(iV + 1) -( 2 -i + 4^ 


w(N) + 1 - 


2N 


w(N — 1) = 0, 


( 2 . 1 ) 


where w(N ) = irG n . An interesting fact is that the formal solution to (12.111 is an asymptotic 
solution; cf. Li and Wong m-, see also m- Hence the asymptotic series (11.511 furnishes a formal 
solution of (|2.1H . Therefore, one way to determine the coefficients a k is to substitute (jl.5p into 
m and equalizing the coefficients of the same powers of x = 1/N. We include some details 
as follows. 


ln(l + x) + y^ 


OO L. 

a k x K 


(i + ®)* 


— 2 - x + 


^2a k x k +( 1 


fe=i 


x\ 2 
2 


ln(l - x) + 


OO u 

a k x 


—' (1 — x) k 
=i y ' J 


= 0 . 


Using the Maclaurin series expansions, we have 

oo oo oo oo / s —2 \ 

— y ^ , s x + y ^ ot k x y ^ d k j~\- k x^ = y ^ I y ^ d k , s OL k d(), s J X =0. 

s =3 k=l j =2 s =3 \k =1 / 

Accordingly, coefficients a k are determined by 

ds— 2 ,s&s —2 H” d s — 3 ?s (T s _3 + * * * d\, s QL\ do ? s — 0, s — 3,4, • • • , 
where the coefficients d s — 2 , s = (s — 2) 2 for s = 3,4, • ■ • , 

, (-l) s + l 1 1 

*- = —-JVT + J^Tj) = and 

j _ ((“1) S fc + l) ( k)s-k (fc) s _fc_i | (k) s -k-2 

^kjS — / 7 \ i / 7 \ i “r 


(s — fc)! 


(s — k — 1)! 4(s — fc — 2)! 


( 2 . 2 ) 

(2.3) 

(2.4) 


for k = 1, 2, • • • , s — 3 and s = k + 3, A; + 4, ■ ■ ■. 

Appealing to (I2.2l) - (l2.4p . the hrst few coefficients a k can be evaluated as 


Ol 4 ) 

a - =_ZL_ 

81.92 ’ 
409875 


«2 192 ; 

— 7615 

«6 — 


8257536 ’ 


a 3 128 ’ 

„ _ 2079 

a 7 = 


oq — 


341 


262144 ’ 


122880 ’ 

_ 679901 

— 1006632960’ 

568756771963 


- _ _ J 6210165 _ 31709469 _ 

9 33554432’ 010 17716740096’ U11 1073741824’ U12 281406257233920' 

One readily sees a periodic phenomenon of the signs of the coefficients, which agrees with 
Theorem |T| To give a full proof of the theorem, we may connect the coefficients with those in 
m, and eventually with a certain hyper geometric function. Indeed, re-expanding the formula 
(11.31) in descending powers of N = n + 1 yields the expansion (11.51) . Hence we have 


Oik = 4 


-k 


1 


— 1 - — 

h ( n 


(k-l)W/3j 
^ (j ~ !) ! ( fc — J') ! 


k = 1 , 2 ,--- 


(2.5) 
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cf. p~3l (4.4)], where ,8j vanish for odd integers j. We also note that the coefficients ( 82 k possess 
a generating function, that is, 


u(x ) = ^ P kx2k 

k =0 


X 

2 sin | 


F 



x 

2 sin | 


F 



( 2 . 6 ) 


where po = 1 and ^ ( 2 s-i)W ’ s = 1, 2, • • • are the positive constants defined in [13], Sec. 3.1]. 
It is shown in m that the generating function u solves a second-order differential equation, 
and consequently the hyper geometric function F(^,|;l;t) is brought in. It is worth noting 
that the function also furnishes a generating relation for the Landau constants, namely = 
oG n x n f° r small x; see [15] . Here and hereafter we denote for short the hyper geometric 
function as F(t) = F (4, 4; 1; t). 


Proof of Theorem [T] From (12.51) we have 

= ^ (~iy +1 Ps ( i N 2fc - 2s 
( 2 k — 1)! “ ( 2 k — 2s)! U 


for k = 1,2, • • • 


Here use has been made of the fact that fas-i = 0 for s = 1, 2, 




(-1) 


fc+i 


Combining (12.61) with 


a 2k x 2k _ 


E 

, s=0 


p s x 


2s 


E 


1 




x 

16 


(2.7) 

. From (|2.7|) we further have 
= u(x)cos—. (2-8) 


and applying a quadratic transformation formula, we have 


1 + E 


(-1) 


fe+i 


- ( 2 k — 1)! 


a 2 k x 2k = 


sin : 


F (sin 2 


1 


sin f cos 2 f 


F ( tan 4 


k =1 v ” r 4 

see [1] (15.3.17)]. Each factor on the right-hand side possesses a Maclaurin expansion with 
positive coefficients; see Nemes m pp. 842-843]. Hence we conclude that 

(—l) fc+1 a 2 fc >0 for fc = 1,2, ■ ■ ■ . (2.9) 

Similarly, we may write 


E 

k =0 


(-l ) fe+1 


( 2 k) 


«2fc+l 2 k 4 J ^ 

- x =zi l^P° x 


2s 


E 


,s=0 




X 


-> = u(x) 

16 1 ' w 


Sill : 


X 


( 2 . 10 ) 


Taking (12.61) into account, we can write the right-hand side term as 


4 cos ■ 


F sin 2 - = 


4 cos | cos 2 § 


F (tan 4 — 


which again has a Maclaurin expansion with all positive coefficients. Here we have used the 
formula 


1 


cos t 


E 

k =0 


t}lE2k f 2k 
( 2 * 0 ! 


where F 2 fc are the Euler numbers such that (—1 ) fc F 2 ^ > 0 for k = 0,1, 2, • • • ; see [171 (24.2.6)- 
(24.2.7)]. Accordingly we have 

(—l) fc+1 a 2fc+ i >0 for k = 0,1,2, • • • . (2.11) 

A combination of (12.91) and (12.111) then gives (11.91) . □ 
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3 Proof of Theorem [2] 

To give a rigorous proof of Theorem 0 we introduce 

R,(N) = e,(N + 1) - (2 - ± + jT) e,(N) + (l - T) E ,(jv - 1) (31) 

for l = 0,1, 2 • • • and N = n + 1 = 1, 2,3, • • • , where £i is the remainder term given in (jl.8p . 

Similar to the derivation of (12.21) . substituting (11.81) into dm and again denoting x = 1/JV, we 
see that Ri(N) is an analytic function of x at the origin, with the Maclaurin expansion 

OO l CO oo 

R i( N ) = ~22 d °’ kxk - ^2 akxk ^2 dk J+ kxJ = *22 r i^ xS ^ ( 3 - 2 ) 

k =3 k=l j =2 s=l-\- 3 

where, for s = l + 3, l + 4, ■ ■ ■, and l = 0,1, 2, ■ ■ ■, the coefficients in (13.21) are 

n, s = ~ (di, s ai + di-i )S ai-i H-h di, s ai - d 0 , s ). (3.3) 

To justify Theorem 0 we state a lemma as follows, leaving the proof of it to later sections. 

Lemma 1. For N = n + 1, it holds 

R 2 i(N):=(-l) l+1 R 2 i(N)>0, n = 1,2,3, - - - , Z = 0,1,2, - - - . (3.4) 

Now we prove the theorem, assuming that Lemma 0 holds true. 

Proof of Theorem 0 For fixed 1,1 = 0,1, 2 • • • , first we show that 

e~2 i(N) := (-1 ) l+1 £ 2 i(N) > 0 (3.5) 

for all n = 1, 2, • • • , where N = n + 1. To this aim, we note that 

S 2 ,(N) = (-1 )' +1 e 2i (iV) = ( ~ 1 01? +1 1 1 + 0 (v) } = fer {' + ' 0 (v) } > 0 (3 6> 

for N large enough; cf. (11.51) . (11.81) and (12.111) . Now assume that e 2 i(N) > 0 is not true for 
some N. Then there exists a finite positive M such that 

M = maxjlV = n + 1 : n € N and e 2 i(N) < 0}. 

Thus for the positive integer M, we have e 2 i(M) < 0, while e 2 i(M + 1), i 2 i(M + 2), • • • > 0. 
Denoting b(N) = (l — 277) 2 for simplicity, from (13.11) we have 

£ 2 i(M + 2) = (1 + b(M + 1 ))£2i(M + 1) + b(M + 1 )(-£ 2 i(M)) + R 2 i(M + 1). 

The later terms on the right-hand side are nonnegative (where M + 1 > 2), hence we have 

£-21 (Tf + 2) > (1 + b{M + l))£2i(M + 1) > £ 2 i(M + 1). 

Moreover, from (13.41) we further have 

£ 2 i{M + 3) > (1 + b(M + 2))s2 i{M + 2) + b(M + 2)(— i 2 i(M + 1)) > i 2 i(M + 2). 
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Repeating the process gives 


£2 i{M + k + 1) > E2i(M + fc), k — 1 , 2, ■ ■ ■ . 

By induction we conclude 

£ 2 i{M + 1) < £ 2 i{M + k ) (3-7) 

for k > 2. Recalling that £ 2 i(N ) = O (lV _2i_1 ) for IV —>• 00 ; cf. (|3.6D . letting k —>• 00 in (13.71) 
gives £ 2 i{M + 1) < 0. This contradicts the fact that £ 2 i(M + 1) > 0. Hence (13.51) holds. 

Now from (jl.8D . (12.91) and (13.51) . we have 

(-1)^2 i-i(N) = + (-l)MAO < 0 for 7 = 1,2,3,--- , and n = 1, 2, ■ ■ ■ , (3.8) 

where IV = n + 1. 

A combination of ((3.5(1 and (|3.8p gives (11.10(1 . Thus completes the proof of the theorem. 

□ 


4 Lemma 12 Estimating of the coefficients ak 


To prove Lemma [T] first we estimate the coefficients a*,, or, more precisely, the quantities 
= ^ ( 2 fc_i y 2k and a 2 k+i = -— ^ ( 2 fc)| 2fc+1 for k = 1, 2, • ■ ■ . These are positive constants; cf. 
(12.91) and (12.111) . As a preparation, we give a brief account of the analytic continuation of the 
hypergeometric function. The reader is referred to m Sec. 3.2] for full details. We denote 

<p(z) =T 1 (sin 2 |) = F Q,^;i;sin 2 ^ for Re 2 € (-2vr, 2vr) U (2 tt, 6tt). (4.1) 

Then the piecewise-defined function 


j ip(z), 0<Rez<27r, 

\ cp(z) ± 2 iip(z — 2 i r), 27 t < Res < 47r and ± Imz > 0 


furnishes an analytic continuation of tp(z) in (14.11) from the strip Rez € [0, 2tt) to the cut 
strip 0 < Rez < 47r and z 0 [2-7T,+ 00 ). What is more, we have the connection formula 
(see [I3l (3.17)]) 

2 

v(z) = va(z) - tp(z — 2ir) In (27r — z) (4.3) 

7r 

for 0 < Re 2 < 47t, with va(z ) being analytic in the strip, and the branch of the logarithm being 
chosen as arg(27r — z) € (—7r,7r). We proceed to show that 


Lemma 2. It holds 


1.9621 1 ^ _ 2.2032 1 

2 k - 1 ( 2 n) 2k ~ ° 2k ~ 2 k - 1 (2vr) 2fc ’ 


k = 9,10,11,- ■■ 


(4.4) 


and 


41n(2jfe + 1) + 0.6551 ^ ^ 41n(2fc + 1) + 2.2048 

- vr(2vr) 2 ^ 1 


k = 9,10,11, - - - . 


(4.5) 














z-plane 


r 



O 



r 


Figure 1: The deformed contour T: the oriented curve (see m Fig. 2]). 


Proof. We understand (12.81) as a generating relation for 012k■ Using the Cauchy integral formula, 
and in view of (12. 6 p . we have 

_ 1 f u(z) cos(z/4) If v(z ) dz 

2k 2ni J z 2 k+i 87 U sin(z/4) z 2k ’ 

where initially the integration path T is a loop encircling the origin anti-clockwise, and is then 
deformed to the oriented curve illustrated in Figure [T] see also m Fig. 2], and v(z) is the 
function defined in (14.21) . From (14.31) . paying attention to the symmetric properties of v(z) and 
r, we have 

~ 1 f v{z) dz 1 f lv(z - 2?r) ln(27r - z) dz 

2k iiri Jf v sin(z/4) z 2k 4iri J Tl sin(z/4) z 2k v h 

where r„ is the vertical part Rez = 37 T, and is the remaining right-half part of T, consisting 
of a circular part around z = 27 t, and a pair of horizontal line segments, respectively along the 
upper and lower edges of (27t, 37t), joining the circle with the vertical line; see Figure [D 
First, straightforward calculation gives 


Iv I < 


M v 1 r°° dy _ M v n 1\ _ 3.1153-•• 

4-7T (37r) 2fc J~oo | sin | 27r(37r) 2fc \4’4/ (37r) 2fc 


(4.7) 


where |u(37t + iy)\ < \Z5maXygR |y?(7r + iy)\ < M v = 2.6393 • • • ; cf. [13, p.297], and B Q, is 
the Beta function. 

Now we turn to the dominant part //. It is readily seen that 


1 f 37r ip(x — 27r) dx 1 f 37r dx 1 f 37r f g(x) — 1 j (x — 27r)^ 

n J 2v sin(x/4) x 2k n J 2n x 2k + it J 2n \ x-2ir } x 2k 


(4.8) 


where g{x) = such that g(27r) = 1. One can see that 9 ^ x _} 2 ^ is positive and monotone 

increasing for x € (27T, 37t] since 


g(x) - 1 = f sin(t/4) 
x — 2tt \tcos(t/4) 


( - 1 

\sin(t/4) 



f --- 

\cos(t/4) 



t = x — 2ir, 
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and each right-hand side term in the curly braces is positive and monotone increasing for 
t € (0,7r]; see pjj p.299] for the monotonicity of ' Therefore, we have for x € (27t, 37t], 


g(x) — 1 g( 37t) — 1 1 

0 < - < —— - := M q = — 

X — 27T 7T 7r 


v/2 F 


1 1 1 
—, —; 1; — 
2 2 2 


- 1 


0.2130-•• . 


Substituting it into (14.81) . we have 


II ~ 2 k - 1 (27r) 2fc + (2k — l)(2n) 2k 

with —3 (§) 2fc < 4fc < 2 ^ g . Further substituting (|4.7I) and (14.91) into (14.61) gives 

- 2_1 4 

a2k 2k - 1 (2vr) 2fc (2k - l)(2vr) 2fc 


(4.9) 


(4.10) 


with 


- {3.1153 (2k - 1) + 3} < 4 < 3.1153 (2k - 1) 


(2\ 2k 2 ttM„ 

W + (FTiy 


Hence for k > 9, we obtain the inequalities in (14.41) . 

Now we turn to the inequality (14.51) for the odd terms. From (|2.6I) and (12.101) we have 


1 f u(z)sm(z/4)dz 1 f v(z)dz 
2k+1 2ni J z 2 k +2 87 H J r cos(z/4:) z 2k+1 ' 

where T is the same path illustrated in Figure [TJ Then, in view of the connection formula (|4.3[) . 
we may write 


_ 1 I" v(z) dz 1 f 

a2k+1 47 n J Fv cos I z 2k+1 + J Tl 

1 f %<p(z - 27r)ln(27r - 
4tt* cos | 


va(z) dz 
cos | z 2k+1 
z) dz 

z 2k +1 


Jv T Ja T Jli 


(4.11) 


where the integration paths r„ and Tj are the same as in (14.61) : see Figured) We note that the 
procedure in da Sec.3.3] applies here, with minor modifications. Case by case estimating gives 


I t | < jMlnf I A) 1 J-H5 3 - , 

'■ vl - 1 2tt V 4 ’ 4 /] (37r) 2fe+1 ~ (37r) 2fe+1 ’ 

see (14.71) . Also, picking up the residue at z = 27r yields 

2u^(27t) 16 In 2 

Ja = ( 2vr ) 2*+l = 7r (27T) 2fc + 1 ’ 


(4.12) 


(4.13) 


where ua(27t) = see da (3.21)]. The dominant contribution comes from the last integral 
J[. We follow the steps in da pp.299-301], and eventually obtain 

j = 41n(2fc + 1) - (47 + 41n(27r)) + 5 itk 1 

1 7r(27r) 2fc+1 ’ ^ ' ’ 
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where |<^fc| < 2M ip + \Mf + rjr e k 2 f° r positive integers k with M v see H3J (3.23)], 

and such that 


ip(x — 27r) 1 

• x-2n x-2n 

bill 4 4 


sm \ J l sin \ 



4 

7T 


for x € (27t, 37t], or, t S (0, 7 r] for t = x — 27r. Here use has been made of the fact that both terms 
in the curly braces are monotone increasing positive functions for t € ( 0 , 7 r]; cf. the derivation 
of (14.91) . Now substituting (14.121) . (|4.13l) and ([4.141) into (14.111) yields 


02fc+l — 


41n(2fc + 1) + (16 in 2 — 4y — 41n(27r)) + 6k 
7r(27r) 2fc+1 


(4.15) 


where |<5fc| < 3.11537T (^y k+1 +2M v ,+ : ^Mf+—^- T e k 2 for positive integers k, which is monotone 

decreasing in k. Straightforward calculation from (14.151) yields ([4.5(1 for k > 9. 

Thus, we complete the proof of Lemma [2J □ 


For later use, we need the following corollary: 


Corollary 3. Assume that otk are the positive constants in Lemma [K Then we have 

^ 2fc < —= 50.8 for k = 5, 6 , • • • , 

« 2fc+2 5 


« 2 fc+l 

«2fc+3 


< 43 for k = 0,1,2,- 


(4.16) 

(4.17) 


and 


(2k + i)^^±i < 0.12 and (2k + \)^±1 < 3.7 f or fc = 0,1, 2, ■ ■ ■ . (4.18) 

OL2k+l d2fc+3 


□ 


The results follow accordingly from Lemma [2] and Table [T| To obtain ([4.17D one may have 
to evaluate the ratio ? 2fc+1 up to k = 13, such that f 2k+1 = 38.578, 38.679, 38.762, 38.829, 38.886 
for k = 9,10,11,12,13. 


k 

1 

2 

3 

4 

5 

6 

7 

8 

9 

a 2k 

a 2k+2 

56.305 

60.184 

57.345 

53.150 

49.797 

47.533 

46.044 

45.031 

44.303 

a 2k-l 

« 2 fc+l 

21.333 

30.720 

34.632 

36.358 

37.227 

37.730 

38.055 

38.282 

38.450 

(2k— l)&2k 
“ 2 fc-l 

0.1041 

0.1184 

0.1007 

0.0851 

0.0749 

0.0684 

0.0642 

0.0612 

0.0590 

(2k—l)a 2 k 

« 2 fc+l 

2.2222 

3.6373 

3.4884 

3.0964 

2.7884 

2.5822 

2.4432 

2.3438 

2.2681 


Table 1: The first few ratios. Calculation conducted using Maple, based on (12.2[ )- (12.4j ). 
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5 Proof of Lemma [T] 

Now that we have proved Lemma [21 we turn to the proof of Lemma [0 
Proof of Lemma [I] To prove (13.41) . the idea is as follows: First we show that 

(- 1)^+2 > 0 for j = l + 1,1 + 2, - ■ ■ , 1 = 0,1,2,--- , (5.1) 

and 

(-1 )* +1 (r 2 i, 2 j+ i + \r 2 i, 2 j + 2 ^ > 0 for j = l + 1, l + 2, • • • , l = 0,1,2, • • • . (5.2) 

Then (13.41) follows immediately from (15.11) and (15.21) since x = 1 /N £ (0,1/2] for N > 2, and 

OO 

R-2l{N) = ^2 (-!) Z+1 {f2l,2j+l + r2l,2j+2 x) X 2j+1 

3=1 +1 

> ( _1 ) i+1 ( r 2l,2j+l + ^ 2 ^+ 2 ) X 2 ° +1 

3=1 +1 V ' 

> 0 


for N = n + 1 = 2,3, ■ ■ •, and l = 0,1, 2, • • • . 

The above idea is simple, yet the verification of (15.11) and (15.21) is quite complicated. We 
begin with ([5.11) . First, a combination of (|2.2j) and (13.3() gives 

(— l)^2Z,2i+4 = d-21+2,21+4 |(~l) Z «2Z+2 j + {“G?2Z+l,2Z+4} | (~l) i+1 «2i+l | > 0 

for l = 0,1, 2, ■ ■ ■ . Here use has been made of (|2.9|) . (|2.11j) . and the facts that d, 2 i+ 2 , 2 i+A > 0 
and d 2 i+i, 2/+4 < 0. Hence (15.11) is true for j = 1 + 1. Therefore, we need only to prove (15.11) for 
j = l + 2, l + 3, ■ ■ ■. In view of (13.31) . it suffices to show, by an induction argument, that 


r 


E 

1,3 


l l-l 

(-l) l+1 22 d 2 k, 2 j+ 2 <+ 2 k > 0 and rfy := (-l) l+1 ^ d 2 k+i, 2 j+ 2 <+ 2 k+i > 0 
k =0 k =0 


(5.3) 


for j = l + 2, l + 3, ■ ■ ■ and l = 0,1, 2, • • • , where «o = — 1. 

Proving rfj > 0 for j = l + 2, l + 3, • • • and l = 0, 1,2, ■ ■ ■ : 

Straightforward verification shows that the first inequality in (15. 3 p holds for 1 = 0: We see 
from (13.31) and (12.31) that 


r 0,j = r 0,2j+2 = 


j + 1 2j + l 


+ > 0 for J = 2 ’ 3 ’ 

8j 


Similarly, from (12.2j) and (13.31) we have rfj = a. 2 d 2 , 2 j +2 + ctodo, 2 j+ 2 - Hence 


r E - 
r l ,3 ~ 


192 


2j + 2 + 


2j - 1 


1 1 , 1 

.7+1 ” 2j + 1 + 87. 


> 


5(j + 1) 

96 


3 (j + 1) 


> 0 


for j = 3,4, • • • . Thus the first inequality in (15.31) is true for 1 = 1. 
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Now assume (15.31) for a non-negative integer l, then, replacing l with l + 2, we have 


r f+2,j = r l,j + «2Z+4 


(2 1 + 3)!d2Z+4,2j+2 


^( 2 1 + 1 )^+ 2 , 2 i+ 2 

OC21+A 


> 0 


for j = l + 4, l + 5, l + 6 , • ■ ■. Indeed, if we write 


(21 + l)!d 2 Z+ 2 , 2 j +2 


( 2 j + 2 / + 2 ) ( 2 j)! ( 2 j — 1 )! 

(2j - 2/)! 4(2j — 2Z — 2)! 


: — Ai + Bi , 


Then, noting that for £ > 0 and j — l > 5, in view of (|4.16l) and Table [TJ we have 


02Z+4 

> (2 j -21- 1)(2 j - 2 l)A t + (2 j -21- 3) ( 2 j -21- 2)B t - 61 (A, + B t ) 

> 90^4/ + 56 Bi - 61 (A l + B,) 

> 0 , 


(5.4) 


since Ai > 4Bi by straightforward verification. Alternatively, applying (14.16|) . for l > 4 and 
j — l > 4, we can modify the above inequalities to give 

ojq/j_9 254 

A+i + B l+l — — - (Ai + Bi ) > 56 Ai + 30 Bi -— (Ai + Bi) > 0 

The remaining cases, namely j = l + 4 with l = 0,1, 2, 3, can be justified by direct calculation: 
The values of 


(21 + 3 )!d 2 z+ 4 , 2 Z+io - ^^(2 1 + l)!d 2 Z+ 2 , 2 Z+io = 62.9, 1004.5, 0.66 x 10 6 , 0.33 x 10 9 , 

«2Z+4 

respectively for l = 0,1, 2,3. Summarizing all above, we see the validity of (15 .4 j) . Therefore, the 
first inequality in (|5.3D is true for j = l + 2, l + 3, ■ ■ • and l = 0,1,2, ■ ■ ■ . 

Proving r^- > 0 for j = l + 2, l + 3, ■ ■ ■ and Z = 0,1,2,---: 

The analysis of ris similar to, and simpler than, that of the even terms rf-. First, for 
l = 0 , the sum in (15.31) is empty and thus we understand that j-q ■ = 0 for all j. Also, it is 
readily seen that rf- = d 1 , 27+204 = ^ for j = 3, 4, ■ ■ ■; cf. (12.41) . Hence, the equality for rfj in 
(15.31) also holds for l = 1. 

Now assume that r f- > 0 for a non-negative integer l and j = l + 2, l + 3, ■ ■ ■ . From (13. 3 p 
we may write 


r ?+2,j — r lj — «2Z+3(2^ + 2 )! d2Z+3,2j+2 + CZ2l+l(2l)\ <^+ 1 , 27+2 := H °7 + c +^l,j (6-5) 

with a positive constant c + = a 2 i+ 3 ( 2 l) \ \d 2 i+i, 2 j+ 2 \, and try to prove that 


(21 + 2 )! 

d<2l+3,2j+2 _ 0121+1 
(2l)\ ( hl+\,2]+2 0-21+3 


6 j + 21 + 2 
6 j + 21 


( 2 j -21- 1 )( 2 j - 21) 


O-21+l 

021+3 


> 0 


for j = l + 4, l + 5, • • • . Here the last inequality comes from (14.171) . Therefore, from (15.51) and by 
induction, we have justified the validity of both inequalities in (15.31) for all j > l + 2 and l > 0. 
Accordingly, we have proved (15.11) for all j > l + 1 and l > 0, noting that the only exceptional 
case j = 1 + 1 has been discussed earlier in this section. 
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Proving (EH): 

In what follows we proceed to prove (15.21) . First, taking into account the formulas (12.21) and 
m, we see that r2i t 2i+3 + b r 2i,2i+4 can be represented as a linear combination of 022+1 and 
022 + 2 - More precisely, substituting in the coefficients c4 jS ; see (12.31) and (12.41) . we have 


(-1 y +i 


1 

122,22+3 + ~1"2l,2l+A 


(21 )\ &21+1 


(2/ + l)(12Z + 5) 
8 


022+2 ( 2 / + 2 ) 2 ( 2 1 + 1 ) 
a2l+l 2 


which is positive for all l since ~ 2 | +2 < ^ 2 i+i) for l > 0; cf. (|4.18l) . Thus (15.21) is true for 
j = l + 1, allowing us to just prove (15.21) for j = l + 2, l + 3, • • • and l = 0,1, 2, • • • . 

For l = 0, it is readily verified from (13.31) and (12.31) that 


1 

_r o,2j+i - 2 r o,2j+2 


1 1 \ / 1 1 

2 j 2j + 2j + \2(2j + l) 4(2j - 1) 



> 0 


for j = 2,3, • • • . Here the right-hand side is the sum of positive numbers when j > 3, and equals 
to ^ when j = 2. Hence (15.21) holds for l = 0. 

For l = 1, recalling that r 2)S = — W 2 ^ 2 ,s + \d\ ts + do,s! cf. (13.31) . from (12.31) and ()2.4I) we may 
write 


1 = _J_\ _ (_± _, JL 

r 2,2j+i + 2 r 2,2j+2 768 + 1536 ^ 2 j 2j + 2 ) \2(2j + 1 ) 4 ( 2 j - 1 )) 16 j ' 


Using the facts that - ^2 - F? for 3 - 3 > and 2 ( 27 + 1 ) _ 4 ( 27 - 1 ) < W, for 3 - ^ we have 
r 2,2j+i + \r 2 , 2 j +2 > — ifj > 0 for all j = 3,4, 5, • • • . Hence (15.21) holds for l = 1. 

Now assume (15.21) for a non-negative integer l, then, from (13.31) we have 


(— l ) l+3 ^r 2 l+ 4 , 2 j+l + -jl'2l++2j+2 S j — (“l ) i+1 ^ r 2Z,27+l + -^21,23+2^ + O t + 
It suffices to show that Oi + Ei > 0 for j = l + 4, l + 5, ■ ■ ■, where for l = 0,1,2,- 

( d 22+3,27+l + 7:^22+3,27+2 ) + 021+1 ( G?22+l,27+1 + 77^22+1,: 


El. 


and 


Oi ■= (-1)' 


Ei ■= (-1)' 


022+3 


«22+4 ( ^22+4,27+1 + 2 ^22+4,27+2 j + «22+2 (<^2+2,27+1 + 7^22+2,27+2 


(5.6) 


(5.7) 


We may write 


Ol = O22+3 + ^2+l) 


«22+l (A + Bi) = O 2/+3 


H;+i + -E>Z+l 


^(Mz+Hz) 

«22+3 


(5.8) 


with Ai = ( 2 J 4 ( 27 ! 1 5 2 Z)! 7;) and Bi = 4 ( 2 ( 7 - 22 2 - 2 )! ■ Observing that A i+1 > ( 2 j - 21) (2j -21- 1 )A t > 
90H/ and Bi+\ = (2 j — 21 — 2) (2 j — 21 — 2>)Bi > 56 Bi for j > l + 5 and l > 0, and recalling that 
? 2i+1 < 43 for l > 0, we have 

cy.o i iq — 


Ol > d 2 2+3 




+ Bi+i 



47 _ ~ 

> 0^022+3^+!, 


j > l + 5, l > 0. 
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Now we turn to £). Similar to the discussion of Oi, we may also write 


Ei = a2i+2(Ci — D{) — a2i+4{Ci+i — A+i), 


where Ci 

1 (2.7—2)! 

4 (2j—2Z—3)! 


_ (2,jf + 1)! _ (2.7-1)! _ (4l+2)(2j—0(2.7 —1)! , f , _ 1 (2.7)! _ 1 (2.7-1)! _ 

{2j-2l)\ (2j—2l—2)\ (2j—2l)\ > dllu 2 (2j-2Z-l)! 8 (2j-2l-2)\ 

. It is readily verified that both constants are positive, and such that |Di < Di < ^Di 


for j > l + 4 and l > 0 with Di 


m 

(2j—2l—l)\ ■ 


Therefore, we have for j > l + 5 and l > 0 that 


47 _ ~ 

Ol + El > —Ol2l+2,Ai + i 

yu 


-Q.2I+2D1 — CH21+aCi + i 


( 2 j - l)!d2; + 3 0 
(2j - 2Z - 1)! 


(5.9) 


where 


fi = (2j - 2/ - 1) 

47 


iL(5i + 7i + 7) - 2(2, f 3fel+4 (2j - l - 1) 
360 «2Z+3 


(21 + l )« 2;+2 j 

«2Z+3 2/ + 1 


> 9 


360 


(5j + 71 + 7) - 0.24(2j -1-1) 


3.7 


-J 


193 . 2077, 2077 
600 J 200 200 ’ 

and thus is positive for j > / + 5 and l > 1. Here use has been made of (14.18p . For the special 
case l = 0 and j > 5, taking Table CD into account, again we have the positivity of H: 


n > (2 j -1) 


47 

360 


(5 j + 7) - 0.24(2j - 1) 


701 

- 2.23 j = -+ 

J 100 1800 


6049.. r , 311 . r , 5 

- 5 > + 900 ^ - 5) 


What remains is the case when j = l + 4 with l = 0,1, 2, • • • . Still we have (15.811 . Since 
A[ + i = i "‘ J 1 4^| 2;+27 ^ > 56H/ and Bi + \ = = 30 Bi, in view of (14.171) we have 


Oi > a.21+3 




13 

56 


13 

(2Z + 7)(12Z + 27) 


Q^Z+S^Z+lj 


from which we see that Oi > ^d2Z+3^.«+i for j = 1 + 4 with l > 2. As a results, we have a 
modified version of (I5.9|) as j = l + 4, 


Oi + Ei > -a.2i+3Ai + i — - 0121 + 2 D 1 — d2z+4C'z+i := —- ~ J+3 H4, 

5 2 b! 


where 


^ _ 12/+ 27 (21 + 3)d2z+4 _|_ 14) (2/ + l)a:2Z+2 Z + 4 


20 «2Z+3 

From (I4.18P we readily see that 


«2Z+3 7(2/+ 1) 


> lAtl! - 0. 12 ( 2I + 14 ) - 3.7 x 1 = 41 + | (I _ 3), 
and is positive for l > 3. 

We £11 the last gap by calculating from (I5.6I) - (I5.7I) that Oi + E[ = 3.3236, 1.9908, 4.3827, 
respectively for l = 0,1,2, with j = Z + 4. Thus we complete the proof of (15.21) . and hence of 
Lemma [D □ 
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6 Discussion 


We have proved the conjecture of Granath |9j, as stated in Theorem [2] and Corollary [TJ of 
which the results of Zhao |20| . Mortici [T3] and Granath [9] are special cases. The asymptotic 
expansion involved, namely (11.51) . corresponds to the special case of Nemes’ expansion (11.21) in 
descending powers of n + h, with h = 1. 

Earlier in na, Li et al consider the case h = 3/4; cf. and (|1.4|) . According to a result 
in m, the error due to truncation is bounded in absolute value by, and of the same sign as, 
the first neglected term for all nonnegative n. As an application, we obtain optimal upper and 
lower bounds up to all orders, holding for all integers n > 0. 

Then, a natural question may arise: 

Question 1. (Li, Liu, Xu and Zhao) Considering the general expansion in (11.21) . for what 
h do we have the “best” approximation in the sense of JlS. Theorem 1] (or, (11.4j) in the present 
paper), or in the sense of Theorem [H and Corollary [7]? 

It is worth noting that the coefficients of the expansion (11.21) possess a symmetric property, 
namely, gk(h) = (—l) fc <7fc(3/2 — h). Hence, if we take h = 1/2, write N = n + 1/2, and specify 
(11.21) as 

OO 

7rG n ~ ln(161V) +7 + X]^- t 6 - 1 ) 

fc=l 1 


Then it is readily seen that 7& = (—1 


, (Z + l)(Z+2) 

and hence (—1) 2 7; +1 < 0 for nonnegative 


integers l, very similar to the result in Theorem [T] Naturally, analysis similar to what we have 

(i+l)(Z+2) 

conducted in the present paper might lead to (—1) 2 e/(lV) < 0 for all n = 0, 1,2, • ■ ■ and 

l = 0,1,2, • ■ ■ , where ei(N) is the error of (16.11) due to truncation at the Z-th order term, with 
N = n + 1/2. 
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